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ABSTRACT 


The  outputs  of  certain  nonlinear  control  systems  are  known 
to  exhibit  jumps  in  their  values  when  the  amplitude  or  frequency  of  the  sinu¬ 
soidal  input  signals  is  varied.  A  semi-graphical  method  is  presented  to  de¬ 
termine  this  jump.  This  method  is  applicable  to  the  class  of  control  systems 
where  there  is  a  single-valued,  piecewise  continuous  nonlinearity  in  the  feed¬ 
back  loop.  The  application  of  this  method  to  a  system  with  a  limiter  is  dis¬ 
cussed.  The  predicted  values  are  verified  by  analog  computer  studies.  Agree¬ 
ment  between  calculated  and  experimental  results  for  third  order  systems  is  very 
good;  for  second  order  systems  the  agreement  is  fair.  Use  of  proportional  deri¬ 
vative  compensation  as  a  method  of  eliminating  jumps  in  the  output  is  also 
examined  both  analytically  and  experimentally. 
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CHAPTER  ONE 


INTRODUCTION 


1  . 1  General 

It  is  well  known  that  the  steady  state  output  of  nonlinear 
feedback  systems  when  subjected  to  a  sinusoidal  input  often  exhibits  one 
or  more  of  the  following  phenomena:  limit  cycles,  jump  resonance  and 
sub-harmonics.  Most  of  these  characteristics,  which  are  associated 
uniquely  with  nonlinear  systems,  have  been  extensively  investigated  by, 
among  many  others,  Hayashi^^^  ^  Ogata Lewis^^\  West,  Douce  and 
Livesley^^\  Prince,  Jr. and  Levinson Many  of  these  investigators 
have  shown  particular  interest  in  the  study  of  jump  resonance. 

1  .2  Method  of  Analysis 

Jump  resonance  is  a  phenomenon  in  which  a  sinusoidal 
input  signal  with  a  fixed  frequency  and  amplitude  may  give  rise  to  out¬ 
puts  of  more  than  one  frequency  or  with  multivalued  amplitudes. 

In  his  paper,  Levinson'  considers  sinusoidal  steady  state 
and  uses  a  semigraphical  technique  to  obtain  the  multivalued  output  and 
the  jump  resonance  conditions  in  a  tachometer-stabilized  system  which 
consists  of  a  linear  plant  and  a  nonlinear  element.  In  essence,  he  equates 
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1  .2  Method  of  Analysis  (cont'd) 

the  equivalent  gain  of  the  nonlinear  element,  expressed  as  a  function 
of  the  error  signal  to  an  expression  involving  the  parameters  of  the  linear 
plant,  the  input  amplitude  and  the  error.  The  latter  expression  represents 
a  family  of  curves  with  w  ,  the  radian  frequency  of  the  input  signal,  as  a 
parameter.  By  plotting  both  sides  of  the  above  equation  separately  on 
the  same  graph,  and  in  terms  of  the  error,  he  is  able  to  determine  the 
jumps  in  the  amplitude  of  the  output  by  noting  the  point(s)  of  intersection 
of  the  family  of  curves  and  the  curves  representing  the  equivalent  gain  of 
the  noni inear  element . 

(2) 

Ogata's  analytic  approach'  '  is  more  interesting.  His 
prediction  of  jump  phenomenon  and  subharmonic  oscillations  is  based  on 
a  method  in  which  the  generally  accepted  Describing  Function  (D.F.) 
gives  way  to  the  Laplace  Transform  technique.  The  Laplace  Transform 
approach  has  a  significant  advantage  over  the  describing  function  method 
in  that  the  assumption  of  a  sinusoidal  input  signal  to  the  nonlinear  element 
which  is  necessary  in  the  describing  function  method  is  no  longer  required, 
and  it  is  therefore  applicable  to  all  types  of  periodic  waveforms  such  as 
square  wave,  triangular  wave,  etc.  However,  the  usefulness  of  this  method 
depends  upon  whether  or  not  the  waveform  of  the  output  of  the  nonlinear 
element  in  question  is  known.  For  instance,  Ogata  assumes  that  the  wave¬ 
form  of  the  output  of  the  nonlinear  element  is  trapezoidal.  The  approxi¬ 
mation  of  the  output  waveform  by  a  trapezoid  may  or  may  not  be  a  good 
one  in  all  cases.  In  view  of  this  basic  limitation,  and  the  fact  that  the 
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1  .2  Method  of  Analysis  (cont'd) 

method  involves  rather  lengthy  and  complicated  mathematical  manipu¬ 
lations, even  in  the  case  of  a  second  order  system,  its  acceptance  over  the 
describing  function  in  solving  nonlinear  control  problems  has  been  limited. 

Another  method,  similar  to  Levinson's,  is  that  due  to  West 
et.al.^'^^  In  this  method  the  authors  employ  the  concept  of  Dual  Input 
Describing  Function  (D.I.D.F.)  and  determine  by  a  semigraph ica I  pro¬ 
cedure  the  occurrence  of  jump  phenomena.  The  region  of  jump  is  indi¬ 
cated  in  the  complex  plane.  The  details  of  this  approach  will  not  be  dis¬ 
cussed  here  due  to  the  limitations  of  space.  It  should  be  mentioned,  how¬ 
ever,  that  this  method  is  even  more  complicated  than  Levinson's  method. 

More  recently  Flatanaka^^^  has  proposed  a  technique  for 
determining  the  frequency  response  of  closed-loop  feedback  systems.  By 
this  method,  the  conditions  governing  jump  resonance  can  be  obtained 
analytically.  Like  almost  all  other  investigators,.  Hatanaka  has  treated 
the  class  of  feedback  systems  in  Fig.  1  .1  in  which  the  nonlinear  element 
is  in  the  forward  path  and  in  series  with  the  linear  part  of  the  system.  The 
closed  loop  system  has  unity  feedback.  Hatanaka  assumes  that  the  non¬ 
linear  element  is  single-valued  and  its  Describing  Function  is  independent 
of  frequency  of  the  input. 


Fig.  1  .1  System  with  nonlinearity  in  forward  loop 
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1  .2  Method  of  Analysis  (cont'd) 

Some  of  Hatanaka's  main  results  are  summarized  in 
Appendix  A.  We  observe  that  E,  the  amplitude  of  the  error  signal  e(t) 
is  the  independent  variable  in  all  the  expressions.  Hatanaka's  approach 
is  not  directly  applicable  to  the  class  of  nonlinear  systems  in  which  the 
nonlinear  element  is  located  in  the  feedback  path  (Fig.  1  .2)  for  two 
reasons.  Firstly,  the  independent  variable  now  becomes  A^,  the  output 
amplitude.  Secondly, the  phase  shift  between  r(t)  and  c^(t)  must  be  taken 
into  account. 


1  .3  Objective  of  This  Thesis 

Since  Hatanaka's  method  is  completely  analytical,  one 
is  likely  to  overestimate  its  usefulness.  However,  practical  implemen¬ 
tation  of  this  procedure  is  difficult  in  view  of  the  complexity  of  the  ex¬ 
pressions  derived.  It  is  the  objective  of  this  thesis  to  present  a  method 
that  is  partly  analytical  and  partly  graphical  for  determining  the  jump 
resonance  frequencies  of  the  class  of  nonlinear  systems  where  the  non¬ 
linear  element  is  in  the  feedback  loop  as  in  Fig.  1  .2.  The  investigation 
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1  .3  Objective  of  This  Thesis  (cont'd) 

wi  II  be  confined  to  systems  in  which  nonlinearity  is  single-valued  and 
whose  describirsg  function  is  independent  of  the  frequency  of  inputs.  The 
nonlinearity  will  be  represented  by  its  Describing  Function  since  we  will 
be  considering  sinusoidal  inputs  to  the  system.  Furthermore,  Describing 
Functions  of  commonly  encountered  nonlinearities  have  been  derived 
and  are  available  in  the  literature  '  '  .  It  is  assumed,  however,  that 

the  llinear  plant  G  acts  as  a  low  pass  filter  and  tends  to  attenuate  the 
higher  order  harmonics  in  the  output  of  the  nonlinear  element.  This 
assumption  is  necessary  (and  hopefully  realized)  to  make  the  Describing 
Function  approach  valid. 
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CHAPTER  TWO 


DERIVATION  OF  CONDITIONS  FOR  JUMP  RESONANCE 


Inpuf-Output  Relationship 


Fig.  2,1  System  with  nonlinearity  in  feedback  loop 

As  mentioned  in  the  previous  chapter,  the  jump  resonance  con¬ 
ditions  will  be  derived  for  a  system  in  which  the  nonlinear  element  is 
in  the  feedback  path.  The  configuration  shown  in  Fig.  1  .2  is  redrawn 
in  Fig,  2.1  for  easy  reference.  G(s)  is  the  transfer  function  of  the  linear 
portion  of  the  plant;  N  is  the  describing  function  of  the  nonlinear  element. 
We  will  restrict  N  to  satisfy  the  following  conditions: 

(a)  N  is  single-valued  and  frequency  independent.  This 
excludes  energy  storage  elements  which  have  hysteresis. 

(b)  N  is  dependent  only  on  the  signal  amplitude  A^  at  its 
input. 

(c)  N  is  twice-differentiable  with  respect  to  A^  .  Conditions 
(a)  and  (b)  are  easily  met  by  many  commonly-encountered 
nonlinear  elements.  Since  most  nonlinearities  in  practice 
are  piecewise  continuous,  or  at  least  their  shapes  can  be 
represented  by  piecewise  continuous  functions,  require¬ 
ment  (c)  poses  no  real  difficulty  in  most  cases. 
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2.1 


Input-Output  Relationship  (cont°cl) 

We  will  determine  the  steady  state  response  c(t)  of  such  a 
system  when  subjected  to  a  single  sinusoidal  input  r(t)o  Let 

r(t)  =  A.  sin  wt  (2.01) 

where  A.  is  the  amplitude  and  w  is  the  angular  frequency. 

In  steady  state,  the  output  of  the  system  will  be 
periodic.  If  the  linear  plant  is  assumed  to  be  a  good  low  pass  filter, 
this  output  can  be  considered  to  be  sinusoidal  and  can  be  represented 
by  its  fundamental  component  or  first  harmonic.  Let  this  component 
be 


c(t)  =  A^  sin  (o)  t  +  0)  (2.02) 

where  A^  is  the  amplitude  and  0the  phase  shift.  This  c(t)  is  the 
input  to  the  nonlinear  element. 

Under  conditions  (a)  and  (b),  the  nonlinear  element 
is  represented  by  its  describing  function  N(A^).  If  c^(t)  is  the  funda¬ 
mental  component  of  the  output  of  the  nonlinear  element,  and 

c^(t)  =  sin  (  (^t  +  0^  (2,03) 


then 


B 

N(A  )  -  -2_ 

°  A 

o 

In  the  frequency  domain,  the  output-input  amplitude 
ratio  of  the  system  is  given  by 


(2.04) 


M  = 


G(a)) 

1  +  G(  oj)  N  (A^) 


A 

o 

Ai 


(2.05) 
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Input-Output  Relationship  (cont'd) 

Let  G(  0))  be  separated  into  its  real  and  imaginary 

parts,  that  is 

G(a))  =  X  ((jo)+  j  Y  (w) 

Substituting  equation  (2.06)  into  (2.05),  we  get 


(2.06) 


=  YOl 


A.  =  ,  v  +  NX)^+  (NY)^  ^ 


V 

Here  N(A^),  X(co),  Y(a)),  are' written  as  N .  X,  Y,  respectively 
for  brevity. 


(2.07) 


2.2  Conditions  for  Jump  Resonance 

With  the  signal  frequency  held  constant,  we  now 
determine  the  input  amplitude  when  output  displays  an  abrupt  jump; 
that  is,  we  want  the  points  at  which 
dA 


8A. 


(2.08a) 


0)  =  const 


or  equivalently, 
9A. 


3A 


=  0 


(2.08b) 
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2.2  Conditions  for  Jump  Resonance  (cont*d) 

These  conditions  are  illustrated  graphically  in  Fig.  2.2. 


A. 


Fig.  2.2  Input-output  curve  with  a  jump 
Using  equation  2.07,  we  find  that 


9  A. 
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0)=  const 


{i/ 


(1  +  NX)^  +  (NY)^ 


-1/2  r 


+  A^  ITo  +  NX)^  +  (NY)^  1  '"•  |(l  +  NX)X-^+  (NY)Y 


dA 


By  setting  this  expression  equal  to  zero,  we  get 

(1  f  NX)^  +  (NY)^  +  A  (X  +  NX^  +  NY^)  =  0 

°  dA 


(2 


or 


X^  (N^  +  Nn'a  )  +  Y^  (N^  +  NN'A  )  +  X(2N  +  N'A  )  +1=0  (2 


where 


dN 


=  N' 


dA 


Equation  (2.11)  gives  the  condition  for  jump  resonance.  This  result  can 
be  further  rearranged  into  a  more  convenient  form: 
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2.2  Conditions  for  Jump  Resonance  (cont'd) 

2 


X  + 


2N  +  N  A 


2(N  +  NN  A^)j 


+  = 


NA 


_2(N^+  NN  AJ. 


This  is  the  equation  of  a  family  of  circles  in  the  X-Y  plane  with 
N  and  A^  as  parameters.  But,  as  noted  earlier,  N  is  a  function 
of  A^,  and  hence  A^  can  be  considered  to  be  the  only  parameter 
in  equation  2,12. 

The  centres  of  these  circles  are  on  the  negative 
real  axis  and  their  coordinates  are 

2N  +  n'a 


2  (N^+  NN  A  ) 
'  o 


0 


The  distance  from  origin  to  the  centre  is 

2N  +  n'a 

CR  =  - 2 — 


2  (N^+  NN  A  ) 
'  o' 


The  radius  of  each  circle  is 


RA  = 


N  A 


2  (N"^+  NN  A  ) 


As  A  takes  on  different  values,  the  radius  and  the  relative 
o  ' 

position  of  the  centres  change.  Fig.  2.3  shows  one  possible 


situation  for  three  values  of  A  . 
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Fig.  2.3  Circles  on  G  plane 

Jump  resonance  conditions  are  indicated  if  the  plot  of  G  (co  )  in  the 
X-Y  plane,  and  the  family  of  circles  in  Fig.  2.3,  intersect.  This  is 


Fig.  2.4  Occurrence  of  Jump  resonance 
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2.3 


Envelope  of  the  Constant  Contours 

For  a  given  G  (o)  ),  it  is  easily  seen  that  there  is 
a  region  in  the  complex  X-Y  plane  where  the  plot  of  G  (  to)  will 
not  intersect  any  of  the  constant  contours.  The  region  in  the 
X-Y  plane  can  therefore  be  divided  into  two  subregions.  If  G  (oo  ) 
lies  in  one  of  these  subregions,  jump  resonance  occurs  and  if  it  lies 
in  the  other,  no  jump  resonance  will  occur.  The  curve  which  separates 
these  two  regions  is  that  which  just  touches  the  constant  A^  circles. 

This  curve  is  the  envelope  of  the  constant  A^  contours,  and  it  can  be 
determined  as  follows: 

It  is  well  known  that  the  envelope  of  a  family  of 


curves  given  by 

f  (x,  y,  c)  =  0  (2.16) 

where  x,  y  are  variables  and  c  is  a  parameter  is  obtained  by  elimi¬ 
nating  c  in  equation  2. 16  and  by  setting 

=  0  (2.17) 

9c 

and  solving  for  c . 

In  the  problem  under  consideration,  we  have 
(equation  2.11)  the  family  of  circles 

(X^  +  Y^)  (N^+ Nn'a^)+ (2N+ n'a  )  X+ 1  =0  (2.11) 

Differentiating  with  respect  to  A  ,  we  obtain 

o  2 

(X^  + Y^)  [3Nn'  + A^  (Nn"+ n'  )1+  (3n'+n"aJX  =  0  (2.18) 

I  I 

where  N  is  the  second  derivative  of  N  with  respect  to  A  . 
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2.3  Envelope  of  the  Constant  Contours  (cont'd) 

For  the  nonlinear  elements  under  consideration, 

i  I  0 

N,  N  and  N  are  all  functions  of  A  ,  and  therefore.  A  can  be 
'  o'  '  o 

eliminated  in  the  equations  (2.11)  and  (2.18).  However,  in  practice 
these  computations  may  be  very  cumbersome  and  a  digital  computer 
may  be  used  to  advantage.  This  will  be  done  in  this  thesis  also. 


2.4 


Input  Amplitude  at  Jump  Resonance  Point 

In  this  section  we  will  use  the  envelope  obtained 

earlier  and  determine  frequency  and  amplitude  of  the  input  signal 

at  which  jump  resonance  occurs.  Figure  2.5  shows  a  case  where  the 

envelope  EN  and  G(  o))  are  shown  intersecting  each  other  at  point  P. 

Let  the  co-ordinates  of  P  be  (X.,  Y.)  and  let  the 

j  I 

associated  frequency  be  oj..  It  must  be  remembered  that  w  is  apara- 


Fig.  2.5  Construction  for  the  location  of  the  centre  of 

a  constant  A  circle 
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2.4  Input  Amplitude  at  Jump  Resonance  Point  (cont'd) 

The  envelope  EN  can  be  considered  to  represent  a 
function  of  X  and  Y.  A  tangent  to  EN  has  a  slope  m  vsrhere 

3Y 


m  — 


3X 


(2.19) 


At  the  point  P,  m.  =  (  —  \ 

«  V  3xA 


(2.20) 


It  is  shown  in  Appendix  B  that  the  equation  of  a 
straight  line  passing  through  P  and  perpendicular  to  the  tangent  at  P  is 


Y  =  -  —  X-tY. +  —  X. 


m. 


I 


I 


(2.21) 


I  I 

This  straight  line  cuts  the  X  axis  at  the  point  X.|  where 

X,  =m.  Y.  +  X. 

■  III 

Let  the  distance  from  P  to  X^  be  L,  then 

L  =  V(V^V^^^  =V<'":  Yj)^+  Y.^ 

Vl 


(2.22) 


I  I 


L  = 


Y. 


(2.23) 


In  equation  2.15,  it  was  shown  that  the  radius  of 
circle  which  represents  the  constant  A^  contour  is  given  by 


RA  = 


N  A 


2(N^+  Nn'a  ) 


(2.24) 


Equations  2.23  and  2.24  can  be  equated  and  solved  for  A^  ,  the 


amplituiie  of  output  at  jump  resonance.  A.  ,  the  amplitude  of  input 

'i 

at  jump  resonance  can  now  be  computed  from  equation  2.07. 
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2.4  Input  Amplitude  at  Jump  Resonance  Point  (cont'd) 

•\/(  1  +  NX.)^  +  (NY.)^ 


A.  = 

I  . 

I 


■v/  X.^  +  Y.^ 

V  I  I 


(2 


2,5  Graphical  Approach 

The  analytical  method,  developed  in  the  previous 
Section,  is  exact.  However,  the  computations  required  in  some  of 
the  steps  are  very  tedious.  For  example,  it  is  quite  difficult  to  ob¬ 
tain  an  analytic  equation  for  the  envelope  expressed  explicitly  in 
terms  of  X  and  Y.  Even  for  a  relatively  simple  nonlinear  element 
such  as  the  limiter,  the  elimination  of  from  equations  (2.11)  and 

I  1 1 

(2.17)  is  difficult  because  N,  N  and  N  are  nonlinear  functions  of 
A^.  In  the  case  of  a  nonlinear  element  represented  by  a  more  compli¬ 
cated  describing  function,  the  analysis  will  be  prohibitively  involved. 
Fortunately,  a  semi-graphical  method  can  be  used.  The  result  ob¬ 
tained  by  this  method  will  probably  be  less  accurate  than  that  obtained 
by  the  analytical  method,  but  the  graphical  method  has  the  advantage 
of  being  simple  and  easy  to  apply. 

We  will  now  outline  this  semi-graphical  method  by 
which  the  input  amplitude  at  jump  resonance  can  be  found  quite 
readily.  The  method  consists  of  five  broad  steps.  Figure  2.6  will  be 
used  to  describe  these  steps. 
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2.5  Graphical  Approach  (cont'd) 

(1)  Both  EN  and  G  (co)  are  plotted  (generally  this  will 
involve  the  use  of  a  digital  computer)  on  the  same 
graph,  and  their  points  of  intersection  P  and  Q  de¬ 
signated  by  (Xp,  Yp,  Wp)  and  (X.2,  Y.2,  a)j2) 
respectively  are  found.  That  part  of  G  ( oj)  lying 
in  the  region  between  P  and  Q  and  bounded  by  the 
envelope  is  the  region  of  jump  resonance.  Operation 
in  this  region  will  cause  a  discontinuous  jump  at  the 

and  tu-o  are 

|2 

the  lower  and  upper  limits  of  the  frequency  of  the  in¬ 
put  signal  between  which  jump  resonance  occur. 

(2)  Within  the  region  of  jump  resonance,  let  D  be  an 
arbitrary  point  where  jump  resonance  occurs  for 
certain  conditions.  Draw  a  tangent  to  G  (co)  at 
this  point. 

(3)  The  straight  line  perpendicular  to  the  tangent  is 

constructed  to  cut  the  X  axis  at  point  X  .  This  then 

c 

is  the  centre  of  the  constant  A  circles  which  just 

o 

touches  G  (w)  at  the  point  D.  We  will  refer  to  this 
as  the  critical  circle. 


output.  The  frequencies  co  p 
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2,5  Graphical  Approach  (cont'd) 

(4)  The  distance  from  the  origin  to  X  is  now  known. 

0 

This  can  be  substituted  into  equation  2.14  and  the 
value  of  A  can  be  obtained  by  solving  the  fol- 
lowing  equation: 


(5) 


CR 


2N  +  N  A 


o. 


X 


2  (N^+  NN  A  ) 
'  o. 


where  N  and  N  are  expressed  as  function  of  A^  , 


Once  A  is  known,  the  corresponding  input  ampli- 
tude  can  be  calculated  from  equation  2,25. 


(2.26) 


Fig.  2.6  Graphical  technique 
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CHAPTER  THREE 


JUMP  RESONANCE  PHENOMENON  IN  A  SYSTEM 

WITH  A  LIMITER 


Introduction 

In  this  chapter,  the  semigraphical  method  described 
in  the  preceding  chapter  will  be  applied  to  a  feedback  control 
system  (Figure  3.1)  in  which  the  nonlinear  element  is  a  limiter. 


Fig.  3.1  System  with  a  Limiter 

Output 

B  ^  Linear  Approximation 


Fig.  3.2  Transfer  characteristics  of  a  saturation  type 
nonlinearity  and  its  I  inear  approximation 
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3.1 


Introduction  (cont“d) 

The  limiter  has  been  chosen  because  it  is  very 
commonly  encountered  in  practice.  Amplifiers  are,  for  example, 
linear  only  over  a  limited  range  of  operation;  a  two-phase  motor 
has  a  maximum  shaft  speed;  and  a  hydraulic  transmission  used  as  a 
motor  device  has  a  maximum  output  speed. 

Figure  3.2  shows  the  transfer  characteristic  of  a 
saturation  type  nonlinearity.  The  limiter  can  be  considered  as  a 
linear  approximation  of  such  a  nonlinearity. 

For  small  signal  inputs  in  the  range  oa  in  Figure 
3.2,  the  device  is  linear  and  output  is  proportional  to  input.  For 
sufficiently  large  inputs,  output  remains  constant  and  becomes  in¬ 
dependent  of  the  magnitude  of  the  input. 

In  Appendix  C,  the  describing  function  (D.F.) 

(3.01) 

It  should  be  noted  that  the  D.F.  is  frequency  in¬ 
dependent  and  as  such  the  procedure  developed  in  Chapter  Two  is 
applicable . 
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Differentiating  N(A  )  in  equation  3.01  with  respect 


to  A^  successively  twice,  we  get 
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3.1 


Introduction  (cont“d) 


or 


and 


N  (A  )  .  A  =  - 
o  o 


B  I 


N  (A^)  =  - 


IT 


'  ■ 


An  examination  of  the  above  expressions  reveals 


(3.03) 


(3.04) 


that  N  (A^)  and  N  (A^)  become  imaginary  when 


-2—  >  1  (3.05) 

A 

o 

It  will  be  recalled  that  equations  (2.11)  and  (2.18) 

I  I  I 

which  determine  the  condition  for  jump  resonance  contain  N  and  N  . 
Consequently,  for  jump  resonance  to  take  place,  the  inequality  3.05 
must  be  satisfied.  This  is  in  accordance  with  the  physical  situation 
at  hand.  When  A^^a,  the  system  is  operating  in  a  linear  region 
with  a  variable  gain  S  in  the  feedback  loop.  In  this  linear  range  of 
operation,  no  jumps  are  possible .  Unless  otherwise  stated,  inequality 
(3.05)  is  assumed  to  hold  throughout  this  thesis. 


3.2  Constant  A  Circles  in  the  X-Y  Plane 

o 

In  this  section  we  will  examine  the  location  of  the 
constant  A^  circles  in  the  X-Y  plane  for  the  system  shown  in  Figure 
3,1  .  We  will  also  examine  the  envelope  of  these  circles. 
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3.2  Constant  Circles  in  the  X-Y  Plane  (cont'd) 


From  Chapter  Two,  we  have 
2N  +  n'a 

CR=  - 2 - 

2N  (N  +  n'a  ) 
o 

where  CR  is  the  distance  of  the  centre  of  the  constant  A 

o 

circles  from  the  origin  in  the  X-Y  plane.  Substituting 
equations  (3.10)  and  (3.03)  into  (2.14)  for  N  and  N  A^, 


we  get 
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Expanding  the  denominator  in  a  power  series,  we  get 
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Clearly,  the  right  hand  side  of  equation  (3.07)  is  always  positive. 

This  indicates  CR,  the  centre  of  constant  A  is  confined  to  the  left 

o 

half  of  the  X-Y  plane. 

Similarly,  the  radius  of  the  constant  A^  circles  is 
obtained  from  equation  (2.15)  as 
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3.2  0>nstant  Circles  in  tfie  X-Y  Plane  (cont'd) 

We  tfien  determine  tf»e  difference  bef^^een  CR  and  RA, 

~2 


CR  -  RA  = 


This  can  be  simplified  into 


CR  -  RA  = 


2B 

.-l/aN 

1 

CN 

o 

"a 

(x-y  V 

'  o  of 

o  ' 

Again,  expanding  the  denominator  in  a  power  series,  (4.10) 
becomes 


(3.09) 


(3.10) 


CR  -  RA  = 


1 


2B 

r  1  /  \  3 

a  1  A  a  A 

1 - 

CO 

o 

o 

-  a 

X"  X  ^X”/  •••  lx-/'  TVA-J 

o  o 

o  o 

(3.11) 


a 

■TT 

o 


Since  ^  1,  we  can  retain  only  first  and  third  onder  terms 


without  much  loss  of  accuracy: 

1 


CR  -  RA 


2B 


2a 

X 


1 


(3.12) 


o  o 

The  right  hand  side  is  always  positive,  therefore 


CR  >  RA 
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3.2  Constant  Circles  in  the  X-Y  Plane  (cont'd) 

This  assures  that  all  the  circles  will  lie  on  the  left  half  of  the 


X-Y  plane . 


If  — <1,  we  can  neglect  the  third  order  term 
A 

o 


also  in  equation  (3.12)  and  we  get 

TT  A 


CR  -  RA  ^ 


4B 


A 

7T  a  o  ^  a 
— . - - - ^ — 

4  B  a  B 


Now  consider  the  limiting  case  when  - . 

A 

o 

The  circle  diminishes  in  size  and  approaches  a  point  since  RA  In 

equation  (3.08)  approaches  zero.  The  difference  between  CR  and 

RA,  in  the  limit,  can  be  obtained  exactly  from  equation  (3.10)  by 

substituting  — ^  =  1  .  Thus, 

A 

o 


CR  -  RA  = 


B 


(3.13) 


(3.14) 


or 


CR  -  RA  = 


since  the  slope  S 


S 

B 

a 


It  can  be  concluded  that  in  the  case  of  a  system 

with  the  limiter  as  the  nonlinear  element  as  shown  in  Figure  3.1, 

all  constant  A  circles  are  confined  to  the  left  of  the  limiting  line 
o  ° 


X=  - 


a 


B 


(3.15) 


(3.16) 
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3.2  Constant  Circles  in  the  X-Y  Plane  (cont’d) 


Fig,  3.3  Limiting  line  and  constant  circles 


The  location  of  this  limiting  line  on  the  X-Y 
plane,  as  shown  in  Figure  3.3,  depends  on  the  slope  of  the 
limiter.  In  the  case  of  a  limiter  with  a  steep  slope,  the  limi¬ 
ting  line  will  be  close  to  the  origin.  It  is  expected  that  inter¬ 
section  (if  any)  with  the  G(a))  locus  will  be  at  the  high  frequency 
portion  of  the  G(w)  curve.  On  the  other  hand,  in  the  case  of  a 
limiter  with  a  small  slope,  the  limiting  line  will  be  quite  far  away 
from  the  origin,  and  intersection  of  G(w)  and  the  envelope  will 
probably  be  confined  to  the  lower  frequency  portion  of  the  G{^) 


curve . 
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3.3  Determ inaHon  of  the  Region  of  Jump  Resonance 

The  semi-graphical  approach  described  in  Chapter 
Two  is  used  in  this  section  to  find  the  region  of  jump  resonance 
in  the  complex  plane.  The  envelope  of  the  constant  circles 
which  forms  the  boundary  of  this  region  is  plotted  with  the  help 
of  a  digital  computer  for  different  values  of  the  slopes  of  the 
limiter.  The  envelopes  for  four  values  of  S  are  shown  in  Figure 
3.4.  It  is  worth  noting  that,  for  values  of  S  equal  to  and  less 
than  1/2,  the  envelopes  are  comparatively  more  linear.  Further, 
the  region  bounded  by  the  envelope  becomes  larger  as  the  slope  of 
the  limiter  increases,  so  that  the  possibility  of  its  intersection  with 
G(a))  also  increases  even  though  the  forward  gain  K  associated  with 
G((jj)  may  be  small.  These  and  other  aspects  of  jump  resonance  will 
be  illustrated  by  means  of  examples  in  the  sections  that  follow. 
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Fig.  3  .4  Envelopes  of  Constonf  Circles  for  Different  Values  of  S 


S=3 


/V 


S  =  5 


Fig.  3.5  Envelopes  of  Constant  A  Circles  and  G  Loci 

o 
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3.4 


A  Second  Order  System 

We  will  choose  a  simple  second  order  plant 
whose  transfer  function  is  G(s),  where 

"  s  (s+  I) 

where  K  is  the  constant  gain  and  is  assumed  to  take  on  different 
values . 

G(to)  =  '  "n  X'"*  ^ — 

(1  +  i  w ) 

In  Figure  3.5,  we  have  drawn  only  the  relevant 
portion  of  the  envelopes  together  with  two  curves  of  G(<^)  when 
K  takes  on  two  different  values.  When  K  =  4,  G(‘^)  intersects 
the  envelopes  corresponding  to  slopes  1,  3  and  5  of  the  limiter. 
When  K  =  2,  G(‘^)  intersects  only  two  envelopes  corresponding 
to  slopes  3  and  5  of  the  limiter.  Therefore,  a  limiter  with  a  slope 
of  1  will  not  cause  a  discontinuous  output  at  any  frequency  if 
K  =  2. 

A  portion  of  Figure  3.5  is  redrawn  in  Figure  3.6 
for  better  clarity,  end  we  shall  make  use  of  it  to  examine  some 
aspect  of  jump  resonance  at  K  =  2. 
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Figo  3o6  Region  of  lump  resonance  for  second  order  system »  K  =  2 
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3.4  A  Second  Order  System^  (cont'^d) 

Region  of  jump  resonance 

By  inspection  of  Figure  3.6,  it  is  seen  that  the 
range  of  frequency  at  which  jump  resonance  will  occur 
is  as  fol lows: 

For  limiter  slope  =  5,  the  range  of  frequency  is 

0.15  (3.19) 

For  limiter  slope  =  3,  the  frequency  range  is 

0.24  ^‘^.^2.2  (3.20) 

Input  amplitude  at  jump  resonance 

We  will  pick  one  point  on  the  G((jj)  curve  lying 
within  the  region  of  jump  resonance  and  determine  the 
input  amplitude  using  the  procedure  outlined  in  section 
2.5. 

At  the  point  P,  (see  Figure  3.6) 

X.  =  -1.6,  Y.  =-3.2  andw.  =  0.5  (3.21) 

I  I  I 

A  tangent  is  drawn  at  P  and  the  line  perpendicular  to 
this  tangent  cuts  the  X  axis  at  -1 ,81  .  The  absolute  value 
1,81  is  substituted  for  CR  in  equation  (3.07).  Thus 
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3.4  A  Second  Order  System  (cont'd) 

Input  amplitude  at  jump  resonance  (cont'd) 

As  assumed  previously,  contribution  from  higher 
order  terms  is  negligible,  so  that  equation  (3.22)  can  be 
simplified  further  as 


1.81  = 


TT-  4“ 

o  Vo 


^  (^)  ] 


(3.23) 


or 

1.81  = 


1  + 


(3.24) 


8  B 
TT'  — 


Y 


where  V  —  for  convenience, 
o 


Consider  the  following  cases: 

B 


Case  1 


slope 
B 


=  5 


Substituting  for  in  equation  (3.24),  we  get 
7.68  -  0.167  =  1 

We  are  interested  only  in  the  real  solution  of  this 

equation  and  it  is  easily  found  to  be 

V.  =  0.515 
I 

Let  a  =  1 ,  then,  since  —  =  0.515, 


A  =1.94 

o. 
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3.4  A  Second  Order  $/stem  (cont*d) 


Input  amplitude  at  jump  resonance  (cont*d) 


From  equation  (3.01) 

N  -  [sin''  V  +  V  yi  - 

Tra  ^ 

(3.28) 

Substituting  for  — ^and  V.  in  (3.28)  we  get 

N  =  3.04 

(3.29) 

Now,  the  input  amplitude  can  be  found  from  equation 
(2.25)  where 


Jo  +  NX  f  +  (NY  )  2 

A.  -  1  1 

A  =  5.7 

(3.30) 

i  lx?  +  Y? 

o. 

1 

^  1  1 

Case  2  slope  — ^  =  3 

In  the  same  manner,  we  obtain  the  following  results: 

V  =  0.478  = 

(3.31) 

o 

A  =  2.09 

(3.32) 

o. 

1 

N  =  1.72 

(3.;33) 

and  A.  =  3.32 

(3.34) 

I. 


I 

A  Third  Order  System 


Consider  the  third  order  plant  with  transfer 
function  G(s)  where 
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3.5 


A  Third  Order  System  (cont'd) 

As  before,  the  envelopes  and  G(‘a))  are  plotted  on 

the  same  graph  paper  as  shown  in  Figure  3.7.  In  this  case,  the 

higher  frequency  portion  of  G(w)  lies  in  the  upper  half  plane  so 

that  the  upper  points  of  intersection  with  envelopes  are  also  in 

this  upper  half  plane. 

Frequency  range  for  jump  resonance 

For  limiter  slope  =  5  0. 125  .7 

For  limiter  slope  =  3  0,23  .4 

^  I 

The  point  we  choose  for  consideration  is  Q  where 

X.  =  -  1.25,  Y.  =  -1.92  0).  =  0.4 

I  I  I 

From  this  we  find  CR  =  1 .42 

Omitting  details,  we  find  the  following  results; 

Case  3  slope  =  5 

-  ^  a 

V.  =  0.56 
I 

A  =1.78 
o. 

I 

N  =  3.25 


and 


Case  4 


A.  =  5.3 

I . 


slope  =  3 

V.  =  0.67 
I 

A  =1.49 

o. 

I 

N  =  2.21 
A.  =  2.96 

I. 

I 


(3.36) 

(3.37) 

(3.38) 

(3.39) 


(3.40) 


(3.41) 


0 
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3.5  A  Third  Order  System  (cont'd) 

The  accuracy  of  the  values  as  obtained  in  sections 
3.4  and  3.5  will  be  verified  when  we  compare  them  with  experi¬ 
mental  values  which  we  will  present  in  the  next  chapter. 
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Figo  3»7  Region  of  jump  resonance  for  third  order  systems.  K  -  2 
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CHAPTER  FOUR 

EXPERIMENTAL  VERIFICATIONS 

Introduction 

This  chapter  is  devoted  to  a  discussion  of  the  ex¬ 
perimental  results  which  were  obtained  by  using  a  PACE  231  R  analog 
computer.  The  systems  considered  in  the  preceding  chapter  were 
simulated  on  the  computer  with  a  view  to  verifying  the  accuracy 
of  the  predictions  made  earlier. 

The  schematic  computer  diagrams  for  the  system 
shown  in  Fig,  3,1  when  G  is  a  second  order  linear  plant  and  when 
G  is  a  third  order  linear  plant  are  shown  in  Fig.  4.1  and  4,2  res¬ 
pectively,  The  details  of  the  simulation  will  not  be  discussed  here 

(12) 

as  the  procedure  is  well  known'  . 


sine  wave  generator  second  order  system  under  study 


Fig.  4.1  Analog  computer  simulation  for  a  second  order  system 
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4.1  Introduction  (cont'd) 


sine  wave  generator 


Third  order  system  under  study 


Fig.  4.2  Analog  simulation  for  a  third  order  system 

4.2  Experimental  Procedure 

The  slope  of  the  limiter  is  set  to  a  specific  value. 

The  frequency  of  the  input  signal  is  kept  constant  at  a  predetermined 
value.  The  amplitude  of  the  input  signal  is  varied  and  the  amplitude 
of  the  output  signal  is  recorded.  All  measurements  are  made  under 
steady  state  conditions. 

Two  values  of  slope  of  the  limiter  S  =  3  and  5  =  5 
are  used  and  the  tests  are  conducted  for  second  and  third  order  systems. 
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4.3  Experimental  Results 

Figures  4.3,  4,4A  and  4.4B  show  the  results 
obtained  experimentally  for  the  second  order  system.  The 
experimental  and  predicted  values  of  input  amplitude  at 
jump  resonance  for  this  case  are  compared  in  Tables  4.1 
and  4.2.  The  actual  waveforms  of  the  system  output  are  shown 
in  Figures  4.5  and  4.6  . 


Table  4<.l  Frequency  range  for  jump  resonance  for  second 

order  system 


Slope  of  limiter 
S 

3 

5 


Experimental 
0.2  <w^3.0 
0.2 


Predicted 

0.24<w^2.2 

0.15<:‘^^2.9 


Table  4.2  Input  amplitude  at  jump  resonance  for  second 

order  system 


Slope  of  limiter 

S 

3 

Frequency  of 
input  signal 
(rad/sec) 

0.5 

Input  amplitude  at 
jump  resonance 
Experimental  Predicted 
3.2  3.32 

3 

1,0 

2.8 

2.82 

3 

1  .5 

2.5 

2.13 

5 

0.5 

5.75 

5.7 

5 

1.0 

5.10 

5.18 

5 

1.5 

4.7 

4.24 

5 

2.0 

4.45 

3.38* 

*  The  predicted  value  atw  =  2  is  not  too  reliable  as  this  frequency  cor¬ 
responds  to  a  point  on  the  G(<jj)  curve  very  close  to  the  X-axis,  thus 
rendering  the  semigraphical  construction  difficult  and  inaccurate. 
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4.3  Experimental  Results  (cont'd) 

Similarly  the  results  for  the  third  order  system  are 
shown  in  Figures  4.7  and  4.8  and  Tables  4.3  and  4.4 


Table  4.3  Frequency  range  for  jump  resonance  for  a 

third  order  system 


Slope  of  limiter 
S 

3 

5 


Experimental 
.23<w  <1.4 
.125<c0^1.7 


Predicted 

.25^^  ^1.5 


Table  4.4  Input  amplitude  at  jump  resonance  for  a 

third  order  system 


Slope  of  limiter 

S 

Frequency  of 
input  signal 
(rad/sec) 

Input  amplitude 
jump  resonance 
Experimental 

at 

Predicted 

3 

0.4 

2.8 

2.96 

3 

1.0 

1.50 

1.59 

5 

0.4 

5.0 

5.3 

5 

0.5 

4.1 

4.7 

** 

5 

1.0 

2.3 

3.2 

**  Footnote  *  applies  also  to  this  case  at  ^  -  1  . 
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Fig.  4a4A  Inpuf-oufput  curves  for  second  order  sysfem,  S  =  5 
Frequency  range  H  CO  ^6 
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Fig.  4.4B  Inpuf-oufpuf  curves  for  second  order  system.  5  =  5 
Frequency  range  0.1  ^<^^0o5 
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Fig.  4.5.  Output  waveform  for  the  second  order  system 
G(s)  =  7^^  O-S.  S  =  5 

when  (A)  A.  =4.5,  (B)  A.  =  6 
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FiQo  4.6  Output  waveform  for  the  second  order  system 

i(^  =1,  S=5 

when  (A)  A.  =  4,  (B)  A.  =  6 
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.  4.  Inpuf-output  curves  for  third  order  system.  5  =  3 
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4.4  Comments 

The  agreement  between  experimental  and  predicted 
vqlues  is  generally  good;  it  is  better  for  the  third  order  system  than 
for  the  second  order  system.  Table  4.2  shows  that,  with  regard  to 
imput  amplitude  at  jump  resonance  in  the  second  order  system, 
agreement  between  experimental  and  predicted  values  is  very 
good  at  low  frequencies,  but  the  discrepancy  increases  as  fre¬ 
quency  increases.  A  rather  large  discrepancy  is  also  observed 
at  the  upper  limit  of  the  frequency  range  (Table  4.1)  for  jump 
resonance  in  the  second  order  system. 

The  discrepancy  between  predicted  and  experimental 
values  of  the  frequency  range  for  the  jump  resonance  is  found  to 
occur  at  all  values  of  slope  of  the  limiter  which  wece  tried.  It  is 
felt,  therefore,  that  there  is  no  apparent  correlation  between  this 
discrepancy  and  the  value  of  slope  of  the  limiter. 

Changing  the  location  of  poles  of  the  transfer 
function  G(s)  to  different  positions  along  the  negative  axis  does 
not  seem  to  have  any  effect  on  the  discrepancy  between  predicted 
and  experimental  values  of  the  frequency  range  and  the  input  ampli¬ 
tudes. 

Based  on  the  above  observations,  it  appears  reasonable 
to  state  that  the  most  probable  explanation  for  the  discrepancies  found 
in  the  case  of  the  second  order  system  is  the  imperfect  filtering  of  high 
frequencies  by  the  plant  G.  The  third  order  plant,  on  the  other  hand, 

t 

offers  a  more  effective  filtering  action,  and  thus  provides  a  better  result. 
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CHAPTER  FIVE 

COMPENSATION  ON  JUMP  RESONANCE 


5.1  Introduction 

In  control  system  design,  the  performance  of  a 
closed  loop  system  must  be  set  to  meet  the  given  specifications. 

One  possible  design  procedure  is  to  develop  first  a  model  system 
satisfying  the  general  overall  requirements,  and  then  take  into 
account  the  detailed  refinements  that  should  be  incorporated  into 
this  fundamental  model  system  so  that  all  the  given  requirements 
are  met.  This  design  procedure  employs  the  commonly  used  tech¬ 
nique  of  compensation. 

In  this  chapter  we  shall  examine  the  advantages, 
if  any,  of  using  feedback  compensation  (derivative  type)  in  order 
to  eliminate  the  jumps  which  would  otherwise  occur.  Elimination 
of  jump  resonance  will  make  the  behaviour  of  the  overall  system 
more  predictable. 

5.2  Possible  Methods  of  Eliminating  Jump  Resonance 

Looking  at  the  curves  in  Figures  3.6  and  3.7,  it  is 
apparent  that  jumps  can  be  avoided  by  shifting  the  G  and  EN  curves 
in  such  a  way  that  they  no  longer  intersect  with  each  other.  This  can 
be  achieved  by  the  relocation  of  either 

(a)  the  G  curve,  or 

(b)  the  EN  curve,  or 

(c)  both  G  and  EN  curves 
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5.2  Possible  Methods  of  Eliminating  Jump  Resonance  (cont'd) 

We  immediately  recognize  that  (b)  is  impractical 
because  the  variables  associated  with  the  nonlinearity  N  are  in 
almost  all  cases  fixed  or  specified.  In  other  words,  it  is  unlikely, 
for  example,  that  the  slope  of  the  limiter  in  a  real  system  can  be 
changed.  Further,  we  may  have  no  control  over  the  point  where 
saturation  level  of  the  limiter  occurs .  Case  (c)  is  eliminated  for 
the  same  reasons.  We  shall  then  explore  further  the  remaining 
possibility,  namely,  relocation  of  the  G{i^)  locus. 


5.3  Effect  of  Varying  K,  the  Gain  of  the  Linear  Plant  G(s) 

In  a  number  of  cases  it  will  be  possible  to  vary 
K  in  G(s)  where  G(s)  is  of  the  form 
G(s)  =  K 

Varying  K  will  change  the  size  of  the  G(t*j)  locus  in  the  X-Y 
plane  but  not  its  shape.  If  K  is  increased,  the  magnification 
factor  M  increases  but  at  the  same  time  the  phase  shift  decreases 
and  vice  versa . 

One  way  of  eliminating  jumps  is  to  decrease  K 
to  a  sufficiently  small  value  so  that  the  G  locus  and  EN  curve  do 
not  intersect.  But  this  has  the  disadvantage  that  the  forward  gain 
will  be  reduced  to  a  rather  small  value  before  the  aim  is  achieved. 


For  the  second  order  system,  for  example, 

K 

s  (s +  T) 
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Figo  5,1  Location  of  G  curves  for  different  values  of  K 
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Fig,  5.  Input-output  curves  for  second  order  system.  K  =  0,6 
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5.3  Effect  of  Varying  the  Gain  of  the  Linear  Plant  G(s)  (cont'd) 

K  will  have  to  be  reduced  to  0.6  before  jumps  can  be  eliminated. 
This  is  shown  in  Fig.  5.1,  and  is  confirmed  experimentally.  The 
experimental  results  are  shown  in  Fig.  5.2  where  it  is  clear  that 
jumps  do  not  occur  when  K  =  0.6.  However,  such  a  small  value 
of  K  may  not  meet  system  specifications  such  as  rise  time,  over¬ 
shoot,  steady  state  error,  etc. 


5.4  Output  Derivative  Feedback 

In  this  section  we  shall  examine  another  approach 
to  eliminate  jumps,  namely,  the  use  of  output-derivative  feedback 
compensation  as  shown  in  Fig.  5.3. 


Fig.  5.3  Output-derivative  feedback  compensation 

The  forward  gain  K  of  the  linear  plant  G  is  held 
constant  while  the  feedback  factor  is  considered  variable.  The 
effect  of  this  feedback  can  be  taken  into  account  by  replacing  G 
by  G  where  G  ,  the  effective  linear  plant,  is  related  to  G  by 

— 1  +  19— 
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5.4 


Output  Derivative  Feedback  (cont'd) 

The  overall  system  is  now  shown  in  Fig.  5.4 


c(t) 


Fig.  5.4  The  equivalent  system 

The  semigraphical  procedure  described  in  Chapter 

I 

Two  can  now  be  applied  to  the  system  in  Fig.  5.4.  The  G  locus 
is  plotted  for  different  value  of  and  superimposed  on  the  EN  curve. 

is  varied  till  the  G  locus  does  not  Intersect  the  EN  curve.  This 
procedure  will  be  illustrated  by  means  of  two  examples. 

Example  1  -  Second  order  system 

"  's  I  +  I) 


and 

g'(s)  =  -  (5.5) 

s  ^  +  s(l  +  2K  ) 

I 

The  G  locus  for  =  0.5  and  the  EN  curve  for  S  =  5  are  shown 
in  Fig.  5.5.  This  system  is  simulated  on  the  computer  and  the 
effect  of  on  the  jump  resonance  is  investigated.  These  results 
are  shown  in  Fig.  5.6.  It  is  clear  that  the  jumps  are  eliminated. 
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Fig.  5.5  Effect  of  oufput  derivative  feedback  on  jump  resonance 
(second  order  system) 
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Fig,  5.6  Inpuf-outpuf  curves  for  compensated  second  order  system 
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Fig,  5,7  Relocation  of  G  curves  by  introducing  compensation 
(third  order  system) 
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5.4 


Output  Derivative  Feedback  (cont'd) 

Example  1  (cont’d) 

I 

It  is  worth  mentioning  that  by  increasing  the  value  of  G 
moves  further  away  from  EN  in  Fig.  5.5.  This  indicates  further 
improvement  in  the  overall  system  performance  as  far  as  avoiding 
the  jumps  are  concerned.  Thus,  a  value  of  K^=0.5  will  insure 
a  smooth  system  response  in  this  case. 

Example  2  -  Third  order  system 

2 

"  - s(s+  1)  (s  +  2) - 


(5.6) 


and 


1  o 

G  (s)  =  - p-2 - - =j - 

sjj  +  3s  +  (2  + 

The  G  and  EN  curves  are  shown  in  Fig.  5.7.  It  is  clear  that  the 
G  curve  must  be  relocated  in  such  a  position  that  the  intersections 
in  both  the  second  and  third  quadrants  of  the  X-Y  plane  are  avoided. 
This  is  accomplished  by  choosing  values  of  such  that  K^“  1  .5. 

It  is  confirmed  experimentally  that  the  input-output  curves  are  smooth 
as  in  the  case  of  the  compensated  second  order  system.  These  curves 
hence  will  not  be  included  here. 
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SUMMARY  AND  CONCLUSION 

A  method  has  been  presented  to  investigate  the 
phenomenon  of  jump  resonance  in  the  class  of  nonlinear  systems  with  a 
single-valued  nonlinear  element  in  the  feedback  loop.  The  Describing 
Function  of  the  nonlinear  element  is  assumed  to  be  frequency  independent.  . 
The  work  presented  here  should  complement  Hatanaka's  work  which  is  suitable 
for  analyzing  similar  phenomenon  in  the  class  of  systems  in  which  the  same 
type  of  nonlinear  element  is  in  the  forward  loop.  The  location  of  the  satu¬ 
ration  type  nonlinearity  will  no  longer,  as  it  has  done  so  far,  cause  any 
difficulties. 

It  appears  that  the  method  presented  in  this  thesis 
is  more  suitable  to  systems  of  third  order  or  higher.  When  it  is  applied  to 
a  second  order  system,  greater  caution  should  be  exercised,  particularly 
in  determining  the  upper  limit  of  the  frequency  range  for  jump  resonance. 

The  reasons  for  this  limitation  have  been  explained  in  the  thesis. 

It  also  seems  possible  that  similar  procedures  could 
be  developed  for  systems  with  multi-valued  and/or  frequency  dependent 
nonlinearities.  Work  along  this  line  is  difficult  but  challenging  enough 
to  warrant  further  investigation. 


-  8S- 


'(  ^  V?*'  i 

HOIZUJ3HCI3  QUA  Y»AMAMJ2 

^  •  I  ^  ^  • 

-  y  t  *e  , 

«Hi  stogitf»vf»l  d  bdnwftiq  n»od  eoH  boHle»tt  A 

o  flliw  ipdftHnon  ^  tK>1p  s»it  fil  930000281  qmuj  I0  nonomonefki 

QfiidhotoO  *<11  .qco^  p^oodbse)  sdt  n}  tofim'plo  iDonllooo  b9ulov*9iQnl2 

« 

.  .Inobnpqdbol  >(0090091^  pd  d  bomuwo  il  foomalo  lotniloon  oHt  noilooul 
siddlut  doidw  ;Ik>w  e'^j^ondoH  Insmalqmoo  bluodt  9i»d  bdri929iq  jHow  offT 
9mo2  odi  doiHw  ni  imoteya  ^0  eooto  pdl  oJ  non^mooodq  loMmii  qnix\iono  lo’f 
-do*  9dl  ncpitoool  9rlT  .qooi  bit>w>cd  odt  nl  21  Inomol*  lopnilfion  9q\1 

4S 

yno  •cuo?  ,io1  oe  onob  lofl  IV  20  ^Toqooi  00  iliw  x^io9fi11non  9q\1  noiloi 

.tdtiuoi^Vib 

i{«9dt  lidt  frf  bdnsediq  boHfom  adi  torb  2iD9qqo  tl  '  ; 

ot  bdiqqo  21  fi  09HW  .ledgM  10  lobio  bildl  ^o  2ffi©^2\2  d  sldoii  u2  enofn  21 
yfiDluoitioq  319X9  sd  biuori^  noiluoo  idoeiQ  ^rnd»y2  lobno  bno392  o 

♦  v» 

,930000291  qitiui  icd  9Qo0i  yaooupoil  ©fit  }o  limU'ieqqu  91+1  gninlrmdsb  ni 
.?l2©dt  ©rif  ni  bsnioJqxs  need  «voH  ooHotifttil  2?i4f  idl  2002091  odT 

•  f 

biuo3  29^win«>Doiq  loimile  torlt  9|'iic2oq  201992  ocio  1l 

'  9 

inabrftqsb  yovi9up9i>  lo\bor>  bouiov-iflufr!  fMw  ^iod2\2  lo^  bdqolovob  od 
dguons  gniQfioilodo  tixi  lluoi^^ib  21  anil  gnolo  ^loW  .eoiliioonitnon 

.noKoQil29vni  ladtiid  inonow  d 


-  59  - 


BIBLIOGRAPHY 


1.  Hayashi,  C.,  "Forced  Oscillations  in  Nonlinear  Systems", 

Nippon  Printing  and  Publishing  Co.  Ltd.,  Osaka, 

Japan,  1953. 

2.  Ogata,  K.,  "An  Analytic  Method  for  Finding  the  Closed 

Loop  Frequency  Response  of  Nonlinear  Feedback 
Control  Systems",  Trans.  AIEE,  Vol.  76  (1957), 

Pt.  11,  pp.  277-285. 

3.  Lewis,  J.B.,  "The  Use  of  Nonlinear  Feedback  to  Improve  the 

Transient  Response  of  a  Servomechanism",  Trans.  AIEE, 
Pt.  11,  Vol.  71  (1952),  pp.  447-453. 

4.  West,  J.C.,  J.L.  Douce  and  Livesley,  R.K.,  "The  DIDF  and 

Its  Use  in  the  Analysis  of  Nonlinear  Feedback  Systems", 
Proc.  lEE,  Vol.  103  (1956),  B,  pp.  463-474. 

5.  Prince,  L.  Jr.,  "A  Generalized  Method  for  Determining  the 

Closed-Loop  Frequency  Response  of  Nonlinear  System", 
Trans.  AIEE,  Vol.  73,  Pt.  11  (1954),  pp.  217-224. 

6.  Levinson,  E.,  "Some  Saturation  Phenomena  in  Servomechanism 

with  Emphasis  on  the  Tachometer  Stabilized  Systems", 
Trans.  AIEE,  Pt.  11,  Vol.  72  (1953),  pp.  1-9. 

7.  Hatanaka,  H.,  "The  Frequency  Responses  and  Jump  Resonance 

Phenomenon  of  Nonlinear  Feedback  Control  Systems", 
Trans.  ASME,  Basic  Eng.  Series  D,  1963,  pp. 236-242. 

8.  Sridhar,  R.,  "A  General  Method  for  Deriving  the  Describing 

Functions  for  a  Certain  Class  of  Nonlinearities", 

IRE  trans,  on  A.C.  Vol.  5,  pp.  135-141,  1960. 

9.  Gibson,  J.,  "Nonlinear  Automatic  Control",  McGraw-Hill, 

New  York,  1963. 

10.  Osgood,  W.A.,  "Advanced  Calculus",  Macmillan  Company, 

N.Y.  1952,  pp.  186-187. 

11.  Jackson,  A,,  "Analog  Computation", McGraw-Hill,  New  York, 

1960. 


\  > 


Y  P 


- 


YWASOOlJSia 


iDi>nilnoH  ni  moilolllazO  baoioT*  ^.D  JrboypH 
^0)^020  ^.blJ  kcO  QniHzildu^  bno  c^n{1ni*i9  noqqiM 

.C29f  ^noqoL 


beeoO  5fft  QnibnH  io^  borhsM  sif^iortA  nA”  ^.X  ^oIdqO 
>l:>odb»d^  >o  senooea^l  y^neupn^  cpoJ 

^(K69r)  6\  .loV  ^331A  .*nofiT  ,'*im»l2x2  loilnoD 

.qq  jr 


sHl  evcnqml  ot  >ldodb«»93  io«nilnoM  »<U  eHT*'  ,,8*L  ^ciwdJ 
^33IA  .eno'iT  \‘*m2inDH^9movi»2  o  senoqzaJI  IneienoriT 

.ee^-^i^.qq  ,(s:e9r)  w  .loV  ji  .is 


boo  3010  •HT*'  ^.X.JI  ^yeUoviJ  bno  aouoO  ,J,l  ,,D.L 
,'’2m9'f2yE  d9odb«93.iD«niinokt  lisytonA  »Hl  nl  ®eU  eti 

^qq  ,a  ,((^f)  GOf  JoV  ,331  ,3oi3 

. 

©Ht  9n»ruofn©1©a  lo^  bofH©M  b9silon©n©0  A**  ^.iC  .J  ^•Dnl'!^ 
,"m©l2y2  loonilnoH  3o  Mnoqi©^  yDn©up9t3^oQj'-b9toiD 

.qq  ,(«<?!)  rr  .f\  .loV;33)A,.2noiT 


rminodsemov^©^  ni  oncHnonsfN  noilDiu1o2  ©mo2"  ,.3  ^noinivaJ 
,''«nal2y2  b«sitidcH2  i»4©moflDoT  ©rft  nonsoriqmH  fhiw 
.9-f  .qq  ,(ee9r)  JoV  ,(f  .13  ,3afA  .20011 


•Doonoiftji  qmuL  bno  lotnoqzaJi  ysnoopaaR  oHT”  ,.H  ^o^lonoloH 
,''2me1*yE  ioilpoD  >lDodb993  toonilnokl  3o  non&monorW 
vvv:-6CS.qq  qC^I  ,Q  esii©^  .cn3  oi2oa  ,3M2A  .enoiT 


pnidiiazaO  ©ril  gniviioQ  io1  boHloM  loi©n©0  A"  ,.3  ,iDHbii2 
^"ioHiiosnilnoU  lo  etolD  niolioD  o  icxi  2noitDnu3 
.0d9f  ,r^f-££f  ,qq  ,S  JoV  .D.A  no  .enoit  331 


^IIiH-woiOdM  ,^loeitaoD  ©ilomotuA  loonilnoU"  ,.L  ,no2d?0 

.€b9i  ,>hoY  woM 


^ynoqmoD  noMimooM  ,'*iuIudIoD  boonovbA”  ,.A.W  ^boogaO 

,’^r-68r  .qq  ,se9r  .y.h 


^^“•oY  woM  %lliH-wt>t'D©M,'*<fcti1otuqfnoD  golonA**  ,.A  ,no2>iooL 

.0d9r 


f*?  ^y  \.^d 


■■ 


'  't 


■  •  f  jg 

i  A  *  .4  -k  ,*1 

‘  ■  i 


-  i: 


^  '1 


-li  »  A 


i'-' 


V*  ■ 


.r 


s 


.£ 


.e 


.6 


'/■I 


.9 


‘it 


■  f 


.J 


-60  - 


APPENDIX  "A" 

SUMMARY  OF  SOME  MAIN  RESULTS  BY  HATANAKA 


Referring  to  Fig.  1.1,  we  get 

r(t)  =  A.  sin  00 1 

I 

e(t)  =  Esin  (<^t+o()  (A.. T) 

c(t)  =  A^  sin  (wt  +  0) 
and  G(‘^)  =  XH  +  iYC*^) 

It  can  be  shown  that  the  output-input  amplitude  ratio  M  is  given  by 

M  =  =  N  (E)  /{X(“)  +[n(E)^‘'}^  +{y(“)15^  (A. 2) 

and  the  jump  resonance  point  is 

|x+  [n(E)]''^  (x+ [N(E)  +  En'(E)]]'H+  =0  (A. 3) 

The  expression  for  the  critical  jump  resonance  curve  is  given  by  the 
equation  (A. 4)  with  E  as  a  parameter. 

X  =  -  n'  (E)  N*  (E)  +  N  (E)  N*'  (E) 

n'  (E)  (n"  (E)}^  +  n"  (E)^N  (E)>^ 

Y  =  +  7n'(e)  n*' (e) '■Cn(e)  -  n*(e)> 

n'  (E)  (n”  (E)^  +  n’7e)(n(E))^ 

N*  (E)  =  N  (E)  +  E  n'  (E) 


where 


A"  XWMWSA 


a. A) 
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(C.A) 


(*-.A) 
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.Ik  m 
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«X  _ 


(3)  H  r  4i-  e  M 
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APPENDIX  "B" 


Referring  to  Fig.  2.5,  slope  of  tangent  to  the  point  P  on  the  curve  G((jo) 
is  given  by 


Equation  of  a  line  perpendicular  to  the  tangent  line  whose  slope  is  m. 


will  have  a  slope  of  -  —  and  the  equation  of  this  line  is 

m . 


(B.l) 


Y  =  X  +  B 

m. 


(B.2) 


where  B  is  a  constant. 

At  w  ,  X  =  X.,  Y  =  Y. 

I  I  I 

then  B  =  Y.  +  — 

I  m. 


X. 


(B.3) 


Substituting  into  (B.2),  we  get 
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Y  =  - 


m. 

I 


X  +  Y.  + 


m. 


X. 


(B.4) 
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(S.8) 


<e.8) 


(k.fl) 
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APPENDIX  "C" 

DERIVATION  OF  DESCRIBING  FUNCTION  FOR  THE  LIMITER 


Fig.  C.  A  limiter  representation 

For  a  linear  element,  transfer  function  defines  a 
unique  input-output  relationship.  In  the  case  of  a  nonlinear  element, 
the  Describing  Function  relates  output  to  input  only  in  some  sense. 
Usually  the  describing  function  of  a  nonlinear  element  is  defined  as 
the  ratio  of  the  amplitude  of  the  fundamental  component  obtained  by 
Fourier  analysis  of  the  output  to  the  amplitude  of  the  sinusoidal  in¬ 
put  (ref.  11),  that  is, 

1^  _  amplitude  of  output  fundamental  component 

amplitude  of  sinusoidal  input 

Let  input  to  the  limiter  be 

R  =  A  sin  wt 
o 


(c.i) 


(C.2) 
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Appendix  "C"  (cont'd) 
Then 

N  =  - 


1 

TT  A 


O  -'O 
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7T 


f  (A  sin  &  )  sin  O'  d  9 
'  o  ’ 


where  9  =  wt 

For  the  limiter  in  Fig .  B .  1 ,  it  is  easily  seen  that 
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where 
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and  substituting  for  S  jn  (C.5),  we  get 
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